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Evaporation is a complex and nonlinear phenomenon due to_ the 
interactions of different climatic factors. Therefore, advanced models 
should be used to estimate evaporation. In the present study, the 
Neural Network-Based Group Method of Data Handling was used to 
estimate and simulate the evaporation rate from the pan in the synoptic 
station of Garmsar city located in Semnan province, Iran. For this 
purpose, the daily meteorological data of evaporation, minimum and 
maximum temperature, wind speed, relative humidity, air pressure, 
and sunny hours of the said station during the nine years (2009-2018) 
were used. The percent of data on training, test, number of the used 
layers, and the highest number of neurons were considered as 60%, 
40%, 5%, and 30%, respectively. The studied method's accuracy was 
investigated using the statistical parameter of Root Mean Square Error 
(RMSE), Mean Absolute Error (MAE) and correlation coefficient, 
and. Sensitivity analysis of the input parameters was performed using 
the GMDH-NN model. This study showed that R’, RMSE, and MAE 
values in the test phase were obtained as 0.84, 2.65, and 1.91, 
respectively, in the most optimal state. From the third layer onwards, 
the amount of the best mean squared errors of the Validation data have 
converged to 0.062, and it is not affordable to use more layers for the 
modeling of the evaporation pan in the Garmsar station. The standard 
deviation and mean amounts of the errors are -0.1210 and 
2.552 respectively. The amounts of the best mean squared errors of the 
validation data are presented. It shows that although the layers are 
increased, the amounts of the mean squared errors have not changed 
considerably. (Maximum 0.003). The sensitivity analysis results 
showed that the two input parameters of minimum temperature and 
relative humidity percent have a higher effect on evaporation pan 
modeling than other input parameters. 
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1. Introduction 


Water is one of the most important factors in developing agriculture and industry in arid and 
semi-arid regions. Abnormal use of water resources and shortage of this vital matter through 
different methods have caused severe challenges and water stress [1]. The increasing population 
and the decline of the natural resources on the planet have forced humans to consider different 
ways to save these resources [2]. 


Water is one of the essential human needs, which has led to more requirements for its planning 
than those in the previous years due to time and space limitations and a very limited volume of 
fresh and exploitable [2]. 


One of the most important water resources meeting the agricultural, drinking, and industrial 
needs in the country, especially in arid and semi-arid regions, is the water stored in lakes behind 
dams and constructed storage pools[3,4]. In several hot and arid areas, a large volume of water 
stored behind dams, agricultural pools, and water storage tanks are wasted due to evaporation 
[5]. Evaporation plays an important role in managing water resources in the region, climate 
change, and agriculture [6]. Considering global climate change, researchers have conducted 
several studies on evaporation worldwide and its assessment for identification in the hydrological 
cycle [7,8]. Evaporation is one of the main phenomena of hydrology [9-11] ,Estimating 
evaporation plays an essential role in estimating the water balance of basins, designing and 
managing irrigation systems, and water resources management. One method for estimating 
evaporation is to use evaporation pans, which are used directly in most parts of the world to 
measure evaporation from the water's free surface [12]. 


The parameters that affect the evaporation rate are relative humidity, temperature, wind speed, 
sunny hours, etc. One of the methods for predicting evaporation is the use of soft computing 
methods. One of the advantages of this method is saving time, reducing trial and error [13,14]. 
Intelligent methods in modeling pan evaporation have been studied and approved by several 
researchers [15—20]. One of the models used in estimating evaporation pan is the developed 
neural network model of the Neural Network-Based Group Method of Data Handling (GMDH- 
NN), between input variables (such as temperature, sunny hours, etc.) and output variable 
(evaporation value). Then, some studies on applying intelligent methods in estimating the rate of 
pan evaporation are referred to. 


J. M. Bruton et al. [21] used the artificial neural network to estimate the daily evaporation from 
the evaporation pan in different parts of the world, including Rome, between 1992 and 1996. In 
this method, precipitation, temperature, relative humidity, solar irradiance, and wind speed were 
used as input data. This study showed that the artificial neural network method has lower error 
than Priestley-Taylor linear multiple regression method. The coefficient of determination and 
root of the artificial neural network model's mean square error was 0.71 and 1.1 mm per day, 
respectively. The R-squared (R2) coefficient and Mean Absolute Error (MAE) of the artificial 
neural network were estimated to be 0.71 and 1.1 mm per day. 
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Keskin and Terzi [22] examined the data of a meteorological station near a lake in the west of 
Turkey to determine the daily evaporation of pan using the artificial neural network model. They 
compared the results of the above designed network model with the results of the Penman 
method. The artificial neural network model results showed a higher correlation with pan 
evaporation rate measured with the Penman method. The best structure of the artificial neural 
network model with 4 input data, including temperature, solar irradiance, air pressure, water 
surface temperature, wind speed, and relative humidity, have a low correlation with evaporation 
intensity in the study area. Qasem et al. [18] predicted the evaporation rate of Tabriz in Iran and 
Antalya in Turkey with three SVR models, ANN, and a combination of them with WSVR and 
WANN wavelet conversion. For both stations, the ANN model has had more reasonable results 
than other presented models had. 


Ashrafzadeh et al. [23] compared the prediction of evapotranspiration in the north of Iran with 
the SARIMA time series model of the intelligent model of support vector machine and Based 
Group Method of Data Handling. The results showed that all three models have good efficiency 
in estimating evapotranspiration. Patle et al. [19] compared MLR and ANN models in estimating 
monthly evaporation of pan in two northern regions of India. The results of this study showed 
that the ANN model had better performance than the MLR model. Alsumaiei [24] modeled daily 
evaporation rate with artificial neural networks in Kuwait. The studied station was Kuwait 
International Airport (KIA). The Meteorological input data of the network include mean 
temperature, wind speed, and relative humidity, which were presented as 4 scenarios.The results 
of this study showed that the combined scenario of mean temperature and wind speed as input 
had better performance than other scenarios in estimating daily evaporation. Al-Mukhtar [25] 
predicted the evaporation rate from the pan in Basra, Mosul, and Baghdad in Iraq. Input 
parameters for artificial intelligence models were minimum and maximum temperature, relative 
humidity, and wind speed. Quantile Regression Forests Model had better performance than 
others. Ashrafzadeh et al. [26] by using the models MLP, SVM, and SOMNN, the pan 
evaporation in Bandar Anzali and Astara, two cities in northern Iran were estimated. Based on 
the high humidity of the studied regions, the input parameters of the three models 
included minimum, maximum, and average temperature, minimum, maximum, and 
average relative humidity, rain, wind speed, and sunshine duration. The SOMNN model worked 
more properly. Singh et al. [27] in another research, the two ANN and MLR models were used for 
the estimation of evaporation. The inputs of the neural network were rain, relative humidity 
delayed for one day, minimum and maximum temperature. The efficiency and 
correlation coefficient of the model ANN during calibration and validation were higher 
than MLR, but the amount of RMSE in the MLR model was higher. 


For accurate prediction with neural networks, a lot of data is needed, so in this study, the gmdh 
neural network was used to predict evaporation from the pan. considering that the GMDH-NN 
intelligent model has not been extensively studied in estimating the daily evaporation rate of pan, 
this study tried to investigate the efficiency of the GMDH-NN model to estimate pan evaporation 
in Garmsar located in Semnan province, Iran. In this study, meteorological data of Garmsar 
synoptic station for 10 years (2009- 2018) were used for modeling. Then, the sensitivity analysis 
of input parameters was performed. 
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2. Methods 


2.1. Study area 


Garmsar city is located in the west of Semnan province in Iran. The distance between this city 
and the capital of Iran (Tehran) is 114 km. The minimum longitude of Garmsar is 51 degrees and 
51 minutes, and the minimum northern latitude is 34 degrees and 18 minutes. The height of the 
meteorological station of Garmsar city center is 899.9 m from sea level. 

Data were analyzed from 2009 to 2018. Data included daily minimum and maximum 
temperature, mean temperature, air pressure, relative humidity, sunny hours, and wind speed that 
were received from the main synoptic station of Garmsar city. The total position of the studied 
station is seen in Figure 1. 
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Fig. 1. Geographical position of Garmsar city 
2.2. Parameters and statistical specifications of data 


In this study, the GMDH-NN model's efficiency for predicting pan evaporation was evaluated 
based on the data on the minimum and maximum temperature, average temperature, relative 
humidity, wind speed, sunny hours, and air pressure, all of which are daily parameters. 


Table 1 shows the studied parameters, abbreviation, and statistical specifications of this research. 
Figure 2 shows the histogram of the input and output data. For better performance, the input and 
output data have been normalized using Relation 1 and Table 2. Thus, all data were between 0.1 
and 0.9 and then used to develop the Relation. This method was used as in the studies [28-30]. 


Table 1 
Statistical Specifications of Data 
Parameter Unit Symbol Mean Standard deviation | Minimum | Maximum 
Minimum temperature ¢ Tinin 13.07 9.71 -12.6 35 
Maximum temperature ¢ Tina 26.31 11.08 -1.6 47 
Relative humidity % RH ean 37.14 19.32 4.5 97.625 
Wind speed m/s WS 7.27 3.99 0 35 
Sunshine hours hr n 8.79 3.27 0 13.8 
Air pressure hPa Pa 914.84 6.26 896.88 936.98 
Evaporation mm E 6.69 7.36 0 39.1 


Parameter — Parameter; 
Parameter, ci.q =| (0.8 = //+0.1 (1) 
Parameter. — Parameter... 
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Fig. 2. Histogram of the input and output data 
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Table 2 
Normalization of the considered data. 
Symbol Parameter Normalized value 
Tinin Minimum temperature ne =0. gf nin — (6.87) +0.1 
a 35.41 
Trix Maximum temperature ‘on =0.8 Prox 2.14 +0.1 
iia 38.26 
RH nean Relative humidity RA oan =0.8 FET aap N98 +0.1 
ii 63.07 
-1.71 
Ws Wind speed WS onnal = je +0.1 
8.55 
n Sunshine hours Normal = 0-8 Ladd +0.1 
8.81 
PA Air pressure PA onal = 9.8 eee +0.1 
16.52 
P- 
E Evaporation normal = 0:8 w +0.1 
Bie 


2.3. Class a evaporation pan 


Class A pan is used in synoptic stations of Iran to estimate the evaporation rate. The evaporation 
data were collected using this pan at the synoptic station of Garmsar. Class A pan is one of the most 
known types of standard evaporation pans used to directly measure the evaporation rate. Internal 
diameter, depth, and water depth are 120, 25, and 20 cm, respectively. The pan has been painted with a 
galvanized sheet. The pan is placed on the wooden bases with a height of 15 cm to be protected against 
heat exchanges with the ground with air rotation below it. Figure 3 shows different parts of the 


evaporation pan Class A. 


Fig. 3. Class A standard evaporation pan. 
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2.4. Neural network-based group method of data handling (GMDH-NN) model 


The GMDH algorithm was first introduced by a Ukrainian scientist named Ivakhnenko [31]. 
GMDH neural network is a self-organizing and unilateral network obtained from several layers, 
each composed of several neurons. All neurons have a similar structure, so that they have two 
inputs and one output, and each neuron with six weights and one bias establishes the processing 
operation among the input and output data based on Relation 2. 

Vue =N(Q, x ;) =p* +w x.) +w ok +w yx — +w ix = +w se ae ae 


7 (2) 


2 
In Relation 2, (i=1,2,3,..., N) where N is the number of observations and (K=1,2,3,.. Em ) and 
ae ih ee where m is the number of the previous layer neurons. 
The weights are calculated based on the Minimum Mean Square Error method and then 
substituted inside each neuron as specified and constant values. 


The obvious characteristic of such a network is that the neurons of the previous phase or the 
2 
previous layer produce new neurons©™ obtained from Relation 3. 


Some of the produced neurons are necessarily removed to prevent network divergence, and 
neurons that remain to expand the network may also be removed due to lack of direct or indirect 
communication with the last layer and creating a network convergence form, called passive 
neurons . The criterion for excluding and selecting a set of neurons in a layer is the Mean 
Square Error (MSE) between the real output and output of each neuron. This criterion for jth 


* 


neuron output, 1.e. es i ), is shown as Relation 4. 


2 m(m —1) 


2 (3) 
iad Ba 
DO, =9 YF 
mse , a cs i 
N (4) 
PEN ie . , 
In the above Relation, j e{ eee m is the number of selected neurons in the 


previous layer. The mapping established between the input and output variables with such neural 
networks as the Volterra nonlinear function is shown as Relation 5. 


A m m m m m m 


y =a,+ > a,%; +) a,x x, +> aaa ee, 
i=l 


i=l j=l i=l j=l k= (5) 


The structure considered for neurons will be in a brief form of two quadratic variables of 
Relation 5: 
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A 
= a 2 2 
y=f (Xj .X ,) = Ay Fax, +4,X , +A,X; +A,X +X xX, 


(6) 


Unknown coefficients®' in Equation 6 are obtained with regression technique such that the 


difference between real output y and calculated values » for each pair of the input variables *: 
“iis minimized. Sets of polynomials are obtained using Equation 6, of which unknown 


coefficients are obtained with the least-squares (LS) method. For each function & (each neuron 
produced), the coefficients of equations of each neuron are obtained to minimize its error in 
order to optimally match inputs with all pairs of input-output sets (Relation 7) [32]. 


0, -G 
= tel > min 


m (7) 


E 


In the GMDH algorithm basic method, all binary compounds (neurons) have made of the n input 
variable, and unknown coefficients of all neurons are obtained using the Least Squares Method. 
Therefore, neurons are made in the second layer according to Relation 8, displayed as Set 9. 


Qpeea 
2 2 (8) 


{o. Xin X ig 
(9) 


We use the quadratic form of the function expressed in Equation 6 for each M of the triple row. 


(i =1,2,...,m) 
p.g €(1,2,...,m) 


These equations can be expressed as the matrix (10): 

Aa=Y (10) 
Where A is the unknown coefficient vector of the quadratic Equation shown in Equation (6), i.e., 
B= NOs, Gy, scp (11) 
And 

Y = Voeus a} (12) 


It can be easily shown from values of the input vectors and function form: 
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2 2 
1 x Pp x qd x Pp x qd x Pp x qd 
2 2 
1 Xoq  Xoqg — XapXrg Xap Xdg 
A= 
1 x x, XX, hi Be 
Mq Mq Mp~~’Mq Mp Mq (13) 


The Least Squares method of the multiple regression analysis yields the solution of equations as 
Equation (14): 


a =(A’A)'A'Y (14) 


This Equation creates the coefficients vector of Equation (5) for the whole M triad set [32]. 

In this study, GMGH neural network was used to model the evaporation pan of Garmsar Station 
for 10 years (6 years for modeling the evaporation pan of Garmsar station and 4 years for the 
model validation). The percent of the Train data, the percent of the Validation data, the number of 
layers used, and the highest number of the used neurons are presented in Table 3. Figure 4 shows 
a schematic view of a gmdh -nn algorithm. 


Table 3 
Specifications of the GMDH Neural Network trained for Modeling Evaporation Pan of Garmsar Station 
P_Train P_Validation Max number of nrouns | Max number of layars 
60% 40% 30 5 


2.5. Modeling accuracy assessment criteria 


Statistical indices of the R-squared coefficient (R’)[33], Root Mean Square Error (RMSE) 
[34,35], and Mean Absolute Error (MAE)[36] were calculated to evaluate the accuracy of the 
intelligent models[37]. The values of these indicators are calculated from the following relations: 


2 


R’= Dial 290 =¥) Z i 
M.-P YO ny 
YO; —x,) 
RMSE =\|i=!—_—_ (16) 
N 
1 n 
MAE = Xo. —*0) (17) 


In relations (15) to (17), *‘is the evaporation rate measured per day, ”' is the predicted 


evaporation rate of the same day, Y is the average values of the measured evaporation, and ’ is 
the corresponding mean for the predicted values. 
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Fig. 4. An illustration of modeling process for the GMDH-NN. 


3. Results 


3.1. Results of GMDH-NN modeling 


In this study, the pan evaporation was calculated every day through the nonlinear GMDH model 
in the Garmsar station. R’, RMSE, and MAE evaluation criteria were used to evaluate the 
performance of this model. Table 4 shows the results of GMDH neural network layers in 
modeling the evaporation pan of the Garmsar station. This table presents the number of neurons 
used in each layer of the trained GMDH neural network and the Mean Squared Error of the 
Validation data for the best neuron of each layer in simulating the evaporation pan of the 
Garmsar station. The results of Table 4 show that the Mean Squared Error of the Validation data 
has converged to 0.062 from the third layer onwards, and it is not cost-effective to use more 
layers in modeling the evaporation pan of Garmsar station. According to Table 5, the GMDH-NN 
model was evaluated for two training and testing phases. The R-squared coefficient of the model 
in the two training and testing phases is 0.86 and 0.84, respectively, indicating that the model has 


a good performance. 


Figure 5 shows the time series of the measured and simulated values with the GMDH model. 
The horizontal axis shows the time series (in terms of the month), and the vertical axis shows 
evaporation values (mm). The more simulated values match the measured values, the more 


accuracy and the less error the model has. 
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Table 4 
Number of layers, number of neurons used, and the best Mean Square Error in GMDH neural network 
trained to model evaporation pan of Garmsar station. 


Layar number Number of nrouns used Best validation RMSE 
1 10 0.065 
2 25 0.063 
3 30 0.062 
4 30 0.062 
5 1 0.062 
Table 5 
The value of error measurement parameters and the accuracy of the proposed model based on GMDH- 
NN. 

R? MAE RMSE 

Training 0.86 1.78 2.49 

Testing 0.84 1.91 2.65 


Figure 6 shows the data predicted by the GMDH-NN model based on the measurement data in 
two training and testing phases. The horizontal axis shows the measured evaporation data (mm), 
and the vertical axis shows the simulated evaporation data (mm). The less the dispersion of data 
around the best fitting line, the more correlation and the fewer errors are achieved. As can be 
seen, the correlation between the measured and simulated fitting data in two training and test 
phases with correlation coefficients of 0.8081 and 0.8598 is relatively high. In addition, Figure 7 
shows that the error values obtained from the developed model based on the group method of 
data handling are small and can estimate the daily evaporation values of the pan. In Figure 8, the 
histogram of the error obtained from the modeling is drawn. The mean and standard deviation 
values of errors are -0.1210 and 2.552, respectively. 
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Fig. 7. Errors obtained from the developed relation in the training and test set. 
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Fig. 8. Simulation error histogram with GMDH-NN 


3.2. Sensitivity analysis results 


GMDH neural network training results showed that after sorting neurons based on the mean 
squares value and removing additional neurons based on standard error, 10 neurons were used in 
the first layer (Table 4). The order of inputs used in the first layer neurons is shown in Table 6. 
Given the results of Table 6, input data 1 and 6, which indicate the minimum temperature and the 
mean relative humidity percent, have the lowest error in simulating the value of pan evaporation 
in Garmsar station, followed by two inputs of the minimum and maximum temperatures with the 
lowest error in simulating the evaporation pan of Garmsar station. The results of Table 6, in 
which the values of the best Mean Squared Error of the Validation data have been presented, 
show that the best Mean Squared Error values have not changed significantly with increasing the 
layers (maximum 0.003). Therefore, two input parameters of minimum temperature and relative 
humidity percent significantly impact the evaporation pan's modeling. Therefore, three input 
parameters of minimum temperature, maximum temperature, and relative humidity percent were 
selected as sensitive parameters for simulation of evaporation pan of Semnan station. 


Table 6 
Mean Squared Error values in the arranged neurons of the first layer of GMDH neural network. 
Nroun number Input data number | Input dataname | RMSE of validation data 

! [1.6] Tininr RH mean 0.0651 
2 [1,2] Trniny Tmax 0.0662 
cS) [1,3] Diag 0.0672 
4 [2,4] Tmaxe WS 0.0674 
5 [1,5] Tmniny PA 0.0676 
6 [2,6] Tinaxe RH mean 0.0680 
7 [2,5] Tinaxy PA 0.0681 
8 [2,3] Tmax, 0 0.0682 
9 [1,4] Tinin, WS 0.0682 
10 [5.6] PA, RH mean 0.0856 
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4. Discussion 


This study aimed to investigate the efficiency of the GMDH-NN model in simulating the value 
of pan evaporation for Garmsar station located in Semnan province, Iran.The study interval was 
from 2009 to 2018. The value of the R-squared coefficient during the test period was 
approximately 0.84. The mean squared error values of the first neuron in each layer of the 
GMDH neural network showed that the RMSE value converged to 0.062 mm after the three 
layers and it is not economical to use more layers in modeling the evaporation pan of this station. 
Also, after sorting the mean squared error in the first layer neurons, two input parameters of 
minimum temperature and mean relative humidity had the lowest RMSE values of 0.0651 and 
were selected as two sensitive parameters in simulating the evaporation pan of Garmsar station. 


Comparing the present study results with the study by Asharfzadeh et al. (2020) [23] shows that 
GMDH neural network model has the necessary efficiency in estimating pan evaporation. 
Karbasi (2016) [38] studied the GMDH model to estimate the evaporation of synoptic stations in 
Ahwaz. The results of this study show that the mentioned model can be effective in estimating 
effective evaporation and can model nonlinear behaviors. In the sensitivity analysis, the results of 
the research by Traore et al. (2010) [39] and Nourani and Sayyah Fard (2012) [40] investigated 
the evaporation using the neural network method. In this research, the most effective parameter 
was the temperature, which is in line with the results of this study. 


The functions of the intelligent methods in estimating evaporation in arid climates are similar to 
the results observed in the previous studies, which have investigated the application of intelligent 
methods for modeling evaporation rates under different climates. Piri et al. (2009) [41] were 
among the first to use ANNs to model pan evaporation rates in arid and semi-arid climates. They 
reported satisfactory performance for ANNs used in a research site located in the southeast of 
Iran. Their study reported R2 = 0.93 for an ANN model with an optimal combination of 4 
meteorological inputs. In the present study, the best obtained R? value was 0.84 during the test 
period, as shown in Table 3. This shows that the models based on intelligent methods are 
effective in arid climates. In addition, the use of artificial intelligence models in arid regions has 
the same prediction error reported in the present study considering the high rate of pan 
evaporation. However, in arid climates, the frequency of such high evaporation rates is higher 
than the pan, leading to lower model performance. 


The present study results are also comparable with the results of other artificial intelligence 
methods used in similar climates. Moghaddamnia et al. (2009) [14] used the ANFIS method on a 
research site located in the Southeast. The R’ value was reported 0.91 for the best performance of 
the ANFIS model during the validation period. However, a similar prediction orientation was 
observed for the ANFIS model. Therefore, further research should improve artificial intelligence 
techniques to allow more reliable predictions for the high rate of pan evaporation. A correction- 
deviation method may indicate a suitable approach in this field. In addition, further studies may 
consider other meteorological variables. Abusada (1988) [42] compared class A pan evaporation 
data collected from Kuwait Airport Station from 1962 to 1977 with theoretical calculations of 
evaporation estimation using the Penman method at the same station for the same period. 
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The comparison showed that the Penman method estimated an annual evaporation rate of 2630 
mm, while the annual evaporation rate of the pan measured for the same period was 3540 mm. 
Therefore, the error resulting from the Penman method is 910 mm per year. Although the 
Penman method is one of the best practical methods for estimating evaporation in an arid climate 
[43], it was found that this method had poor performance compared to the artificial neural 
network models. The evaporation pan wall prevents additional sunlight and increases heat 
exchange with the surrounding atmosphere [44].Therefore, physical models cannot be used to 
estimate the evaporation of the pan directly. 


5. Conclusions 


In this study, estimation of pan evaporation values in Garmsar city located in Semnan province, 
Iran, was investigated using the developed neural network model of group method of data 
handling. In this research, R’, RMSE, and MBE evaluation criteria were used to evaluate the 
model results. The results of this study show that the GMDH-NN model is suitable for modeling 
"evaporation pan processes". Also, the mean squared error values in the first neurons of the 
trained GMDH neural network layers showed that the RMSE value of the first neuron in the third 
layer converges to 0.062 mm and it is not necessary to use more layers. The mean squared error 
values in the first layer neurons showed that the two input parameters of minimum temperature 
and relative humidity have the lowest RMSE values (0.0651 mm) in simulating the amount of 
evaporation pan of Garmsar station and were selected as two sensitive parameters. Arid and 
semi-arid climates have unique climate regimes that are characterized by "scarce water 
resources, "bare vegetation," and "high evaporation rate ".Considering the report of the Food and 
Agriculture Organization of the United Nations (FAO), excessively arid climates are defined as 
the regions where annual precipitation does not exceed 3% of the annual evaporation. Comparing 
the performance of the ANN model with other practical models is essential by dealing with the 
subject of contribution of pan wall to heat exchange. Then, it can be noted that considering that 
the time behavior of daily evaporation is non-stationary, it is better to investigate and compare 
other intelligent methods for modeling. 
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